A New Galactic Extinction Map in High Ecliptic 

Latitudes 



Tsunehito Kohyama 1 , Hiroshi Shibai 1 , Misato Fukagawa 1 , Takahiro Sumi 1 , 

and 

Yasunori Hibi 2 

1 Department of Earth and School, Graduate School of Science, Osaka University, 1-1 
Machikaneyama-cho, Toyonaka, Osaka 560-0043, Japan 
kohyama@iral. ess.sci. osaka-u. ac.jp, shibai@iral. ess.sci. osaka-u. ac.jp, 
misato@iral. ess.sci. osaka-u. ac.jp, sumi@iral. ess.sci. osaka-u. ac.jp 
2 Advanced Technology Center, National Astronomical Observatory of Japan, 2-21-1 Osawa, Mitaka, 

Tokyo 181-8588 
yasunori.hibi@nao.ac.jp 

(Received ; accepted ) 
Abstract 

In this study, we derived a galactic extinction map in high ecliptic latitudes for 
| /5 | > 30°. The dust temperature distribution was derived from the intensities at 100 
and 140 fim with a spatial resolution of 5'. The intensity at 140 fim was derived from 
the intensities at 60 and 100 fim of the IRAS data assuming two tight correlations 
between the intensities at 60, 100, and 140 /xm of the COBE/DIRBE data. We 
found that these correlations can be separated into two correlations by the antenna 
temperature of the radio continuum at 41 GHz. 

Because the present study can trace the 5 '-scale spatial variation in the dust tem- 
perature distribution, it has an advantage over the extinction map derived by Schlegel, 
Finkbeiner, and Davis, who used the DIRBE maps to derive dust temperature dis- 
tribution with a spatial resolution of 1°. We estimated the accuracy of our method 
by comparing it with that of Schlegel, Finkbeiner, and Davis. The spatial resolution 
difference was found to be significant. The area in which the significant difference is 
confirmed occupies 28% of the region for | j3 \ > 30°. 

With respect to the estimation of extragalactic reddening, the present study has 
an advantage over the extinction map derived by Dobashi (2011), which was based on 
the 2MASS Point Source Catalog, because our extinction map is derived on the basis 
of far-infrared emission. Dobashi's extinction map exhibits a maximum value that is 
lower than that of our map in the galactic plane and a signal-to-noise ratio that is 
lower than that of our map in high galactic latitudes. This significant difference is 



confirmed in 81% of the region for | /3 \ > 30°. 

In the areas where the significant differences are confirmed, the extinction should 
be estimated using our method rather than the previous methods. 

Key words: ISM: dust, extinction — infrared: ISM — galaxies: ISM — methods: 
data analysis 

1. Introduction 

Radiation from extragalactic objects is reduced by galactic interstellar dust, whose ef- 
fect becomes appreciable in the ultraviolet (UV)-optical wavelength region. The foreground 
emission by the dust is significant in the cosmic background radiation at the infrared-radio 
wavelength (CIB; CMB). Therefore, the galactic dust component must be subtracted along 
each line of sight for the study of the cosmic background radiation. 

For this purpose, Schlegel, Finkbeiner, & Davis (1998; hereafter SFD) created an all- 
sky galactic extinction map by using the all-sky far-infrared maps of the Infrared Astronomy 
Satellite (IRAS) and the Diffuse Infrared Background Experiment (DIRBE) on the Cosmic 
Background Explorer (COBE). The spatial resolutions of the IRAS and DIRBE data are 5' 
and 0°.7, respectively (IRAS Explanatory Supplement; DIRBE Explanatory Supplement). SFD 
derived the dust temperature from the DIRBE data at 100 and 240 /im because the far-infrared 
intensity at wavelengths shorter than 100 fim includes the emissions from very small dust grains 
that are transiently heated to high temperatures because of the absorption of high-energy 
photons (Aannestad & Kenyon 1979; Draine & Anderson 1985; Desert, Boulanger, & Puget 
1990). SFD derived an extinction map using IRAS at 100 fim intensity and dust temperature 
derived from the DIRBE data. Therefore, this map may not trace dust variations in the scale of 
a few minutes, as reported by Arce Sz Goodman (1999). SFD calibrated the dust column density 
value into the color excess by using the observed reddening of elliptical galaxies cataloged in 
Faber et al. (1989). 

Several extinction maps have been derived on the basis of star catalogs for dark cloud 
detection. Dobashi et al. (2005) derived a visual extinction map with a 6' spatial resolution by 
using the Digitized Sky Survey (DSS). Dobashi (2011) derived an all-sky extinction map with 
a V - 12' spatial resolution using the 2 Micron All Sky Survey Point Source Catalog (2MASS 
PSC). In addition, he derived the color excess of a star by using the "X percentile method," 
which is a recent extension of the Near Infrared Color Excess method (NICE, Lada et al. 1994; 
Rowels & Froebrich 2009). Compared with the NICE method, the X percentile method is more 
robust against contamination by foreground stars because it statistically corrects for the effects 
of stars in the foreground. He has renewed the estimation of the background star colors using 
a star catalog simulated by Robin et al. (2003). These methods measure the dust only in front 
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of the stars and are based on near-infrared radiation. Therefore, the precision of these methods 
can degrade in directions of high extinction. In addition, the signal-to-noise ratio (S/N) can 
decrease in directions of small extinction, such as that in high galactic latitudes. Therefore, 
these methods may be ineffective for the correction of galactic extinction in extragalaxies. 

We created a galactic extinction map of the Cygnus region (80° < I < 90°, —4° < b < 8°) 
using a new method (Kohyama et al. 2010) that derived the dust temperature distribution with 
a 5' spatial resolution. We obtained the intensity at 140 /mi from the IRAS data at 60 and 100 
/an, assuming two tight color-color correlations identified between DIRBE data at 60, 100, and 
140 /mi. The spatial resolution of our dust temperature map was 10 times higher than that of 
SFD, and our map could detect the galactic dust behind the star shown in the DSS catalog. 

In this study, we apply our method to a region of high ecliptic latitudes for | /3 \ > 30°. 
In contrast to what was previously reported (Kohyama et al. 2010), a high galactic latitudes 
region is included here. The extinction map derived in this study is compared with those of 
SFD and Dobashi (2011). The difference between this study and SFD in the high galactic 
latitudes region is expected to be smaller than that in the galactic plane because the dust 
distribution in high galactic latitudes is more uniform than that in the galactic plane. The 
difference between this study and Dobashi (2011) is considered to be smaller than that between 
our previous study and Dobashi et al. (2005) because the near-infrared study in Dobashi (2011) 
allows the detection of a larger extinction using the color excess of more distant stars than that 
in the optical study by Dobashi et al. (2005). 

The Far-Infrared Surveyor (FIS) on the AKARI satellite has a spatial resolution less 
than 1' (Murakami et al. 2007; Kawada et al. 2007). FIS has observed the all-sky intensities at 
65, 90, 140, and 160 /im, and its diffuse intensity maps are scheduled for publishing. The High 
Frequency Instrument (HFI) on the Planck satellite has observed the all-sky intensities at 350 
/im to 3 mm (Planck HFI Core Team 2011a, b). Planck collaboration (2011) derived the dust 
temperature with a 4' spatial resolution. A comparison of the independent results of the FIS, 
HFI, and this study is expected to improve the accuracy of dust temperature and extinction 
measurements. 

The far- infrared data and our analysis are shown in Section 2, the color-color correlations 
are determined and separated in Section 3, the maps of dust temperature and extinction are 
derived in Section 4, the result of this study is compared with those of SFD and Dobashi (2011) 
in Section 5, and the conclusion is provided in Section 6. 
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2. Analysis and Data 



2.1. Analysis 

In general, interstellar space in the galaxy is optically thin in the far-infrared wavelength 
region. The far-infrared intensity is written by Equation (1). 



where I U (X) is the intensity at wavelength A, t wo Mm is the optical depth at 100 /im, /3 is the 
spectral emissivity index, B u (X,Td) is the Planck function, and T d is the dust temperature. 
We adopted (3 = 2 according to the previous studies (Draine & Lee 1984; Hibi et al. 2006; 
SFD, Kohyama et al. 2010). In addition, we adopted a one-zone model such that the dust 
temperature is assumed to be constant along each line of sight, according to the previous 
studies (SFD; Kohyama et al. 2010). 

Color-correction should be applied to the observed data with the photometric bands of 
DIRBE and IRAS because their published all-band intensities were derived by assuming that 
the source spectrum ul u (\) is constant. The conversion from the observed (color- uncorrected) 
intensity to the color-corrected intensity is given by Equation (2) with a color-correction factor 



determined from the spectral emissivity index and dust temperature if the spectrum is given 
by Equation (1). The color-correction factor is unique for each instrument, and those of the 
DIRBE and IRAS intensities are referenced in the respective explanatory supplement. In this 
study, the color-corrections are applied for the intensities at 100 and 140 /mi, but not for that 
at 60 /j,m, as reported in previous studies (Hibi et al. 2006; Kohyama et al. 2010). This is 
because the spectral energy distribution shape around 60 /^m cannot be described by Equation 
(1). This correction affects the intensities by a few percentage points, and is therefore not 
crucial for our results. 



The far-infrared data are summarized in Table 1. All data are presented in the 
Hierarchical Equal Area isoLatitude Pixelization (HEALPix) projection (Gorski et al. 2005). 
We set the pixel scale in this study to be as large as 1'.72, which corresponds to Nside = 2048; 
Nside is the pixel scale parameter of the HEALPix projection. The HEALPix pixel scale has 
previously been defined discretely (Gorski et al. 2005, Table 1). The Improved Reprocessing 
of the IRAS Survey (IRIS) by Miville-Deschenes & Lagache (2005) was used as the IRAS im- 
age. The IRIS images are available in the HEALPix projection with a 1'.72 pixel scale. The 




(1) 



of K x (P,T d ): 




, jC (A) represents the col 



2.2. Data 
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Zodi Subtracted Mission Average data (ZSMA) was used as the DIRBE image, as reported in 
Kohyama et al. (2010). The ZSMA images are available in the Quadrilateralized Spherical Cube 
(QSC) projection, and the data are binned to 393,216 pixels with a 19'. 5 pixel scale. The direct 
conversion from the QSC to HEALPix 1'.72 pixel size was found to be inappropriate because it 
produced blank pixels and forced interpolations during the conversion process. Therefore, we 
first attempted to convert the QSC data to the HEALPix projection with a 27'. 5 pixel scale 
(Nside = 128). Conversion to the larger pixel size was not a disadvantage for this study because 
the ZSMA images need to be spatially smoothed to achieve sufficient S/N (see Section 3-1 for 
details). The images were smoothed by Gaussian convolution with a 30' Full Width at Half 
Maximum (FWHM) to reduce the jagged pattern caused by the conversion. The la noise level 
of DIRBE140 is 2.4 MJysr -1 in the QSC projection (DIRBE explanatory supplement), and has 
accordingly become 1.4 MJysr -1 after the processes described above. 

The interplanetary dust (IPD) emission in both IRIS and ZSMA images was first sub- 
tracted on the basis of the model of Kellsal et al. (1998). Then, the CIB offsets at the three 
wavelengths were subtracted in this study. The CIB intensities employed in this study (shown 
in Table 2) are lower than the upper limits indicated by Hauser et al. (1998). 

Table 3 shows the notations for IRAS and DIRBE intensities at 60, 100, and 140 ixm. 
The IRAS intensity at 140 fim was derived in this study. 

3. Color— Color Correlations 

Hibi et al. (2006) analyzed the ZSMA data in the galactic plane (| b | < 5°) and identified 
two color-color correlations at 60, 100, and 140 fim. Kohyama et al. (2010) analyzed the 
ZSMA data in the Cygnus region and identified two correlations that differed slightly from 
those identified by Hibi et al. (2006). In this study, we analyzed the ZSMA data and identified 
two correlations in the region for | j3 \ > 30°. The high galactic latitudes region was included 
in this analysis. The S/N of DIRBE140 is lower than those of DIRBE60 and DIRBE100. In 
addition, the S/N of DIRBE140 in high galactic latitudes is lower than that in the galactic plane. 
Therefore, the DIRBE maps were smoothed spatially by Gaussian convolution to improve the 
S/N of DIRBE140. Two correlations are found, and we separated the two correlations by 
using the Wilkinson Microwave Anisotropy Probe ( WMAP; WMAP Seven- Year Explanatory 
Supplement) all-sky map at 41 GHz produced by Jarosik et al. (2011). 

3.1. Smoothing of DIRBE maps and Identification of Color-Color Correlations 

As shown in Table 4, we changed the FWHM of the Gaussian depending on the S/N of 
DIRBE140 because the maps smoothed with a large beam have the disadvantage of throwing 
away small-scale information in the high-S/N regions. 

Figure 1 shows the S/N distribution of DIRBE140 before the processes described in 
Section 2-2. The FWHM of the Gaussian for the smoothing is defined to achieve an S/N of 
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DIRBE140 over 10 (Table 4), with which S/N Hibi et al. (2006) found tight correlations in 
color-color diagrams. 

Figure 2 shows the color-color diagrams in the region for | f3 \ > 30°. Two correlations 
with the different inclinations are seen in the color-color diagrams. The galactic and magellanic 
cloud data points shown in the figure are color-corrected. We define the regions of the LMC 
and SMC as those for which IRAS100 is over 10 MJysr -1 and 3 MJysr" 1 in the directions of 
LMC and SMC, respectively. 

3.2. Division of Color-Color Correlations 

Figure 3 shows two color-color diagrams with antenna temperatures both lower and 
higher than 0.7 mK at 41 GHz. The former exhibited a slower inclination than the latter. In a 
manner similar to that reported by Kohyama et al. (2010), we described the former and latter 
as slow and steep cases, respectively. The correlation coefficient of the slow and steep cases are 
—0.6 and —0.8, respectively. The two color-color diagrams appeared to show tight relations. 
Hibi et al. (2006) separated the two correlations by using the 10 GHz radio continuum emission 
map (-4° < I < 56°, | b | < 1°.5) produced by Handa et al. (1987). Hibi et al. (2006) and 
Hirashita, Hibi, & Shibai (2007) indicated that the inclination of the correlation in the color- 
color diagram changes depending on the intensity of the free-free emission on the line-of-sight. 
Hibi et al. (2006) divided all data points into the two groups with the border of the 0.5 K antenna 
temperature at 10 GHz. The correlation lower and higher than 0.5 K is called as the main and 
sub correlation, respectively (Hibi et al. 2006). The 10 GHz radio continuum emission map 
traces the free-free emission, but is available only within a small part of the sky. Therefore, 
we employed the 41 GHz radio continuum map for separating the data points into the two 
correlation groups because it also traces the free-free emission that is available for the whole 
sky. Kohyama et al. (2010) separated the two correlations such that DIRBE60/DIRBE100 is 
0.28. A specific value of DIRBE60/DIRBE100 cannot separate the two correlations found in 
this study. 

3.3. Fitting 

Hibi et al. (2006) and Kohyama et al. (2010) fitted the color-color correlations as a 
power-law function represented by 

DIRBEU0 _ / DIRBEQ0 \ b 

DIRBE100 ~ a \DIRBE100J ' ^ ' 

We also fit the color-color correlation in Equation (3) using the least squares method with 
parameters of coefficient and power in the two color-color diagrams (Figure 3); the results are 
summarized in Table 5. We tried two cases for the fitting: for case I, each data is given a 
weight according to its la photometric error, while for case II, no weight is applied. According 
to Kohyama et al. (2010), only the power- law index b will be used, but the coefficient a will 
not be used in the following discussion (Equation 4). 
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The fitting results are shown in the color-color diagrams in Figure 3. As seen in the 
right panel of Figure 3, the slopes of the power-law function for the steep case are similar 
among cases I and II and previous studies. However, the difference between cases I and II is 
significant for the slow CclSGS, £lS shown in the left panel. Case I represents the correlation of 
the DIRBE data with a relatively high S/N (over 10 for DIRBE140). Case I is consistent with 
the main correlation (Hibi et al. 2006) and the slow case derived in Kohyama et al. (2010). The 
two latter cases were derived for the galactic plane. It is reasonable to expect that case I is 
consistent with the two latter cases because high S/N data concentrate to the galactic plane. 
On the other hand, case II weighs on the data with smaller S/N (but it is still large enough), 
which are mostly seen above the galactic plane. We add the difference between the results for 
cases I and II to the uncertainty of our method because it is difficult to accurately determine 
the spatial boundary of cases I and II from the present data. Case II may hold above the 
galactic plane, while case I may hold at the galactic plane. 

3.4- Uncertainty of the Interplanetary Dust Emission Subtraction 

The appropriate subtraction of the IPD emission is difficult to achieve because the IPD 
emission alone has not been observed. To avoid an uncertainty of the IPD emission subtraction, 
we did not analyze the region for | /3 \ < 30°, which shows the systematic stripe features appeared 
in DIRBE60 (and IRAS60) in the ecliptic plane (DIRBE60; Hauser et al. 1998; IRAS60; Miville- 
Deschenes & Lagache 2005). The stripe features were parallel to the ecliptic plane and did not 
appear in the region for | (3 | > 30°. Kellsal et al. (1998) predicted the IPD emission with 1% 
precision; however, there is significant uncertainty regarding the stripe feature. For example, 
the uncertainty is about 50% at the ecliptic longitudes of 180°. The average of DIRBE60 
adjacent to the stripe is 0.6 MJysr -1 and inside the stripe is 0.3 MJysr -1 . 
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4. Galactic Extinction Map 



We derived IRAS140, dust temperature, and extinction for cases I and II. The calculation 
method was the same as that reported in Kohyama et al. (2010), except for the power-law 
values of the color-color correlations. The power-law values were determined on the basis of 
the antenna temperature at 41 GHz of the line of sight. IRAS140 was derived by the following 
equation: 

IRASUO _ DIRBEUO / IRASQO . DIRBEQ0\ b 

IRAS100 ~ DIRBE100 \IRAS100 ~ DIRBEIOOJ ^ 
The dust temperature was derived from intensities at 100 and 140 /mi, and the optical depth 
at 100 /mi was derived from the dust temperature and the intensity at 100 /mi. The optical 
depth was converted to the extinction at the V band (i.e., A v ) assuming the extinction law of 
Mathis (1990). 

7 I/C (140 \im) of IRAS is the color-corrected intensity, which is derived as described below. 
7^(140 /im) of IRAS is derived by the color-uncorrected IRAS and the color-corrected DIRBE 
data by Equation (4). From these 7^(140 /tm) and 7„(100 /im) of IRAS, the dust temperature 
is calculated once. The dust temperature is used for the color-correction of 7j,(100 /im) of IRAS 
to obtain 7^(100 /xm) of IRAS. Finally, 7^ c (140 //m) of IRAS is derived from 7^(100 fim) and 
7„(60 /im) of IRAS by Equation (4), and the dust temperature is derived by 7^(100 \im) and 
7 VjC (140 nm) of IRAS. 

The notations for Ay shown in this paper are summarized in Table 6. Figures 4 and 5 
illustrate the maps of dust temperature and extinction, respectively. 

We employ the results based on case I for the following reasons. The correlations in 
case I are believed to represent a common dust property described below. The color-color 
correlations of case I are similar to the correlations identified by Hibi et al. (2006) and those 
of extragalaxies. Hibi et al. (2006) found that the color-color correlation of the magellanic 
clouds were similar to the main correlation, which is the slow case of Hibi et al. (2006). We also 
found that the slow case I correlation is similar to the main correlation. Moreover, Hirashita 
& Ichikawa (2009) found that the color-color correlations of the blue compact dwarf galaxies 
observed by AKARI are similar to the correlations identified by Hibi et al. (2006). 

The difference between the results for case I and the color-color correlations found by 
the previous studies are not added to the uncertainty of our method. It is because that the 
difference between the results for cases I and II is larger than that between the results for case 
I and Hibi et al. (2006) and that between the results for case I and Kohyama et al. (2010) in 
both slow and steep cases. As demonstrated in Section 5, the Ay difference between SFD and 
the result for case I and that between Dobashi (2011) and the result for case I are significant 
even if the difference between the results for cases I and II is added to the uncertainty of our 
method. A significant difference is also expected even if the difference between the results for 
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case I and Hibi et al. (2006) and that between the results for case I and Kohyama et al. (2010) 
are added to the uncertainty of our method. 

5. Results and Discussion 

As described above, we have derived the maps of dust temperature and extinction with 
a 5' spatial resolution on the basis of far-infrared intensities. The derived extinction map based 
on case I was compared with the maps of SFD and Dobashi (2011). 

5.1. Comparison with SFD 

The difference in the spatial resolution of the dust temperature is expected to create 
dispersion in the difference between this study and SFD, but it should not cause a systematic 
difference. We estimated our accuracy assuming that SFD's map is correct in a 1° or larger 
spatial scale. In addition, we assumed that the difference in spatial resolution of the dust 
temperature causes dispersion only in the difference. Figure 6 shows the comparison between 
the Ay values of case I and SFD, in which dispersion difference is evident. Ay (SFD) is derived 
from the E(B — V) of SFD assuming Ry = 3.1. In particular, a systematic difference appears 
for Ay (SFD) > 1 mag. The mean and median of the Ay (case I) j Ay (SFD) ratio for Ay (SFD) 

< 1 mag are 1.05 and 0.96, respectively. On the other hand, Ay (case I) is systematically larger 
than Ay (SFD) for Ay (SFD) > 1 mag. The mean and median of the Ay (case I) /Ay (SFD) 
ratio for Ay (case I) > 1 mag are 1.21 and 1.13, respectively. From the assumption described 
above, we estimated the accuracy of Ay (case I) to be as much as 5% and 21% for Ay (SFD) 

< 1 mag and Ay (SFD) > 1 mag, respectively. 

The reason why Ay (case I) differs systematically from Ay (SFD) for Ay > 1 mag is 
described below. The region for Ay > 1 mag corresponds to lower galactic latitudes regions 
(| b | < 15°). In these regions, the dust temperature varies in a small spatial scale and can become 
quite low and high in limited portions of the sky. High-temperature regions with scales of a few 
minutes cannot be resolved in the DIRBE maps. In contrast to the DIRBE maps, the IRAS 
maps resolve these high-temperature regions. Therefore, the IRAS dust temperature map shows 
a larger area of the high-temperature region when compared to the DIRBE dust temperature 
map. In these regions, the scatter diagram between the dust temperature distributions derived 
by the DIRBE maps and that derived by the IRAS maps will show that the dust temperature 
derived by the IRAS maps is systematically lower than that derived by the DIRBE maps. 
When the dust temperature decreases, the dust extinction increases. Therefore, the extinction 
derived by the IRAS maps becomes systematically larger than that derived by the DIRBE 
maps. SFD derived the dust temperature using the DIRBE maps, but we derived the dust 
temperature using the IRAS maps. Therefore, in these regions (i.e., Ay > 1 mag), the dust 
extinction derived by our method is systematically larger than that derived by SFD. 

The difference in the spatial resolution of the dust temperature contributes to the hor- 
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izontal dispersion in Figure 6; this trend changes at Ay (SFD) = 0.1 mag. The horizontal 
dispersion for Ay (SFD) < 0.1 mag appears relatively wider than that for Ay (SFD) > 0.1 
mag. The wider dispersion may be caused by SFD's assumption that dust temperature is uni- 
form in the small extinction region (i.e., in high galactic latitudes), and they therefore combined 
the DIRBE maps with background averages obtained at high galactic latitudes for | b | > 75°. 
No such assumption was made in our method. 

Figure 7 shows a histogram of Ay (SFD) vs. Ay (case I). The mean, median, and 
standard deviation are —0.28, —0.001, and 2.29 mag, respectively. Figure 8 shows the Ay 
difference distribution of Ay (SFD) vs. Ay (case I). The absolute Ay difference between Ay 
(case I) and Ay (SFD) in high galactic latitudes is smaller than that in the galactic plane 
because of the smaller extinction in high galactic latitudes. In addition, the dust distribution 
in high galactic latitudes is more uniform than that in the galactic plane, and the standard 
deviation in high galactic latitudes is significantly smaller than that in the galactic plane. The 
standard deviation is 0.10 mag for | b | > 45° and 2.66 mag for | b \ < 45°. 

Figure 9 shows the uncertainty of our method estimated by the following equation: 

AA V = A{noise) + A(case I — case II) + A(accuracy of the case I). (5) 

The first term is the Ay la noise error calculated from the IRIS la noise at 60 and 100 /im. 
Photometric noise could cause the horizontal dispersion in Figure 6 in addition to the difference 
in the spatial resolution of the dust temperature. The second term is the difference between 
Ay (case I) and Ay (case II) calculated as | Ay (case I) — Ay (case II) | . The third term is the 
accuracy of Ay (case I) as estimated by the systematic difference between this study and SFD. 
The values of the third term are 0.05A V (SFD) and 0.21A V (SFD), for Ay (SFD) < 1 and Ay 
(SFD) > 1 mag, respectively. This is a dominant term for Ay (SFD) > 1 mag (| b |< 15°). The 
second term is dominant for Ay (SFD) < 1 mag (| b |> 15°). Near the galactic pole, the first 
term becomes equivalent to the second term because of the low S/N of the IRAS data. The 
extinction value should be estimated by our method rather than by SFD's method when the 
difference from Ay (SFD) is larger than AAy given in Equation (5). Figure 10 shows the area 
in which the Ay difference between this study and SFD is larger than the Ay uncertainty of 
our method, as shown in the following equation: 

| Ay(SFD) — A v (case I) | > AA V . (6) 

The area represented in Figure 10 occupies 28% of the region for | /3 \ > 30°, and exists in 
not only the galactic plane but also high galactic latitudes. The area in which the significant 
difference is confirmed increases if we consider just case I or II because the nominal uncertainty 
in each case is smaller than the difference between Ay (case I) and Ay (case II). If we assume 
case I only, the area occupies 49% of the region for | /3 | > 30°. 
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5.2. Comparison with Dobashi (2011) 

Ay (2MASS) is converted from E(J — H) derived by Dobashi (2011), assuming the 
extinction law of Cardelli, Clayton, & Mathis (1989) as reported in Dobashi (2011). E(J — H) 
is the most sensitive quantity in the color excesses derived by Dobashi (2011). The case of 
E(J — H)x m , used for the dark cloud survey in Dobashi (2011), is compared in this study (see 
the definition of E(J — H)^ m described in Dobashi 2011). The LMC and SMC were masked 
because Dobashi (2011) did not derive the extinction of the LMC and SMC. 

Figure 11 shows the result of the comparison between this study and Dobashi (2011). 
The A v (2MASS) are saturated at Ay (2MASS) > 20 mag, indicating that the method of 
Dobashi (2011) has an upper limit of Ay (2MASS) = 20 mag. It is apparent this method is 
not useful for Ay > 20 mag even though the upper limit of Ay = 20 mag is four times as 
large as that of the visual extinction derived by Dobashi et al. (2005). The visual extinction 
derived by Dobashi et al. (2005) was saturated at Ay > 5 mag in the Cygnus region (Kohyama 
et al. 2010). Ay (2MASS) is systematically smaller than Ay (case I) for 1-10 mag because 
galactic dust may exist behind the stars cataloged in the 2MASS PSC. Dispersion that appears 
in the vertical direction for the Ay (case I) < 1 mag is probably caused by the low S/N of Ay 
(2MASS) in high galactic latitudes. In high latitudes, Dobashi (2011) employed a beam size 
that was larger than that in the galactic plane to achieve a constant noise level. The extinction 
in high latitudes was smaller than that in the galactic plane. The la noise level is 0.4 mag in 
Ay (2MASS). The mean uncertainty of our method was 0.07 mag for | b \ > 45°. 

Figure 12 shows the histogram of Ay (2MASS) — Ay (case I). The mean, median, 
and standard deviation of the Ay difference are —0.58, —0.09, and 3.69 mag, respectively. 
Figure 13 shows the distributions of the Ay difference between this study and Dobashi (2011). 
As expected in the systematic difference in Figure 11, a systematic difference is seen in the 
galactic plane. The mean and median for | b \ < 45° are —0.75 and —0.10 mag, respectively. 
Though galactic dust may exist behind the stars cataloged in the 2MASS PSC, Ay (2MASS) 
is systematically larger than Ay (case I) toward the galactic bulge. It indicates that Dobashi 
(2011) may overestimate the extinction toward the galactic bulge. A scatter appears in high 
galactic latitudes. The standard deviation for | b | > 45° is 0.34 mag. 

Figure 14 shows the area in which the difference in Ay obtained in this study and 
Dobashi (2011) is larger than the Ay uncertainty of our method, as shown in the following 
equation: 

| A V (2MASS) - A v {case I) \ > AA V . (7) 

The area represented in Figure 14 occupies 81% of the region for | j3 \ > 30°. The Ay (2MASS) 
and Ay (case I) are in agreement at middle-galactic latitudes, which indicates that the extinc- 
tion derived at near-infrared wavelengths above the middle-galactic latitudes may trace most 
of the dust column in our galaxy along a line of sight. 
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In the area represented in Figure 14, the extinction should be estimated using our 
method instead of Dobashi (2011)'s method with respect to the estimation of the extragalactic 
reddening. However, our method may not be suitable as an estimator of the reddening for 
objects located in the galaxy because the extinction derived in this study reflects the entire 
integration through the galaxy. Note that the reddening derived by Dobashi (2011) is larger 
than that derived by our method toward the galactic bulge as shown in Figure 13. 

The area in which the significant difference is confirmed increases if we consider just case 
I or case II. If we assume case I only, the area occupies 90% of the region for | ft \ > 30°. 

The Ay comparison results are summarized in Table 7. 

Future prospects are described below. Observation results of AKARI/FIS should be 
compared with the results of this study. FIS has observed the all-sky intensities at 65, 90, 140, 
and 160 fim with a spatial resolution less than 1' (Murakami et al. 2007; Kawada et al. 2007). 
The FIS results are independent of the results of this study. The FIS observed data at 65 and 
100 /im will be calibrated by the IRAS and DIRBE data. The FIS observation data at 140 /im 
can be calibrated by the DIRBE data. Therefore, with respect to the spatial resolution, the FIS 
calibration accuracy at 140 /im is 1°. The accuracy is expected to increase upon comparison 
with this study. Planck/RFI showed the dust emission distribution with a 4' spatial resolution. 
The spatial resolution was nearly equal to that obtained in this study. Planck collaboration 
(2011) derived the dust temperature distribution on the basis of intensities at 100 and 350 
/im. The IRAS data were used as the intensities at 100 /im in the Planck collaboration (2011). 
Planck collaboration (2011) assumed a spectral emissivity index of (5 = 1.8. If a significant 
difference between Planck's result and this study is confirmed, the validity of the spectral 
emissivity index of the dust emission should be examined. We will produce maps of the dust 
temperature assuming a spectral emissivity index of f3 = 1.8. 

6. Conclusion 

We derived a new galactic extinction map in high ecliptic latitudes for | /3 | > 30° on the 
basis of the far-infrared intensity and dust temperature with a spatial resolution of 5'. The dust 
temperature was derived on the basis of intensities at 100 and 140 /im. The intensity at 140 
/im was derived from the intensities of the IRAS data at 60 and 100 /im, assuming color-color 
correlations between the intensities of the DIRBE data at 60, 100, and 140 /im. 

We identified two tight correlations in the region for | /3 | > 30°. The DIRBE maps were 
spatially smoothed to increase the S/N of DIRBE140, particularly in high galactic latitudes. 
We found that the two correlations can be separated by the antenna temperature of the radio 
continuum at 41 GHz. 

Our method had a clear advantage over SFD because the spatial resolution of the dust 
temperature was 10 times higher than that of SFD. The accuracy was estimated to be as much 
as 5% and 21%, for Ay (SFD) < 1 and Ay (SFD) > 1 mag, respectively. A significant difference 

12 



in Ay was confirmed in the comparison between this study and SFD. The area in which the 
deviation was larger than the Ay uncertainty of our method occupied 28% of the region for 
| (3 | > 30°. In this area, the extinction should be estimated using our method instead of SFD's 
method. 

The extinction map derived in this study was compared with that derived by Dobashi 
(2011). The A v (2MASS) was saturated at Ay > 20 mag. A low S/N of Ay (2MASS) 
appeared in high galactic latitudes. The S/N of Ay (2MASS) was significantly smaller than 
that determined by this study in high galactic latitudes. The area in which the significant 
difference was confirmed occupied 81% of the region for | /3 \ > 30°. In this area, with respect to 
the estimation of extragalactic reddening, the extinction should be represented by our method 
as opposed to Dobashi (2011)'s method. 

All analyses described in this paper were carried out in Interactive Data Language (IDL) 
with the HEALPix package (K. M. Gorski et al. 2005, ApJ, 622, p759). The COBE/DIRBE and 
WMAP Seven- Year data were obtained from the Legacy Archive for Microwave Background 
Data Analysis (LAMBDA). The IRIS data were obtained from the IRIS project, and the SFD 
data were obtained from the Dust Maps project. The 2MASS extinction data were obtained 
from the Astro Laboratory of Tokyo Gakugei University. This study was supported by Grant- 
in- Aid for challenging Exploratory Research (22654023). 
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Table 1. Summary of the far-infrared data. 





Wavelength 


resolution 


1(7 

noise level 


Calibration 
error 


1PD emission 
model 


KlYnopf inn mnlhnn 

used in this study 


Pixel scale 
used in this study 




H 


[arcmin] 


[MJysr- 1 ] 


[%] 






[arcmin] 


IRAS 


60 


4'.0 


0.03 


10.5 








IRAS 


100 


4'.3 


0.06 


13.5 


Kellsal et al. 






DIRBE 


GO 




0.09 


10.5 


(1998) 


HEALPix 


1'.72 


DIRBE 


100 


42' 


0.1 


13.5 








DIRBE 


140 




2.4 


10.6 









IRAS data are the IRIS data (Miville-Deschenes, & Lagache 2005). The DIRBE data are ZSMA (DIRBE explana- 
tory supplement). 



Table 2. CIB intensities at 60, 100, and 140 fim. 



Wavelength [/im] 


Intensity [MJysr 1 ] 


Reference 


60 


0.27 ± 0.05 


Miville-Deschenes, Lagache, & Puget (2002) 


100 


0.73 ± 0.03 


Planck collaboration (2011) 


140 


0.94 ± 0.16 


Matsuura ct al. (2011) 



Table 3. Notations for the IRAS and DIRBE data. 

IRAS60 intensity observed by IRAS at 60 /xm 

IRAS100 same as above, but at 100 /im 

IRAS 140 intensity derivied in this study at 140 /im 

DIRBE60 intensity observed by DIRBE at 60 /im. 

DIRBE100 same as above, but at 100 /im. 

DIRBE140 same as above, but at 140 /im. 
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Table 4. Summary of the spatial smoothing process by Gaussian convolution used to achieve the S/N of DIRBE140 over 10. 



S/N ratio 
before tne smoothing 


FWHM of a Spatial resolution 
Gaussian [deg] after the smoothing [deg 


Expected S/N 
ratio improvement 


> 10 




0.9 


X 1 


5-10 


1.5 


1.8 


x 2 


2.5 - 5 


3.4 


3.6 


x 4 


1.3 - 2.5 


7.0 


7.2 


x 8 


< 1.3 


13.8 


14.4 


x 16 




1 ^^^= 5 

Fig. 1. Map of the S/N of DIRBE140, which was divided into five classes. The white region represents 
classes with S/N ratios higher than those of the black region, and the gray zone represents the region that 
is for /3 | < 30°. The distribution is shown in the Mollwcide projection with the galactic coordinate. The 
dotted lines indicate the galactic latitudes of | b \ = 45° 
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DIRBE140 / DIRBE100 DIRBE140 / DIRBEIOO 

Fig. 2. (left) Color-color diagram of the region for | /3 | > 30°. The DIRBE maps were smoothed to 
achieve an S/N of DIRBE140 over 10. The dots and squares represent galactic and magellanic data 
points, respectively, (right) The contour map illustrates the density of the data points on the diagrams. 
The bin size is set to 0.02 in logarithmic units. Contour lines are drawn according to the number of data 
points, including 0.5 (to express the boxes that have only one data point), and are 0.1N, 0.5N, N, 2N, 3N, 
and 5N on the diagram, where N = 40 indicates the average number of data points contained in each bin 
on this contour map. 



18 




DIRBE140 / DIRBE100 DIRBE140 / DIRBEIOO 



Fig. 3. Two color-color diagrams of the region for (3 | > 30°. (left) 41 GHz antenna temperature below 
0.7 mK; (right) 41 GHz antenna temperature over 0.7 mK. Solid and dashed lines show the color-color 
correlations for cases I and II, respectively; dotted lines show the result in the galactic plane (Hibi et 
al. 2006), with the main correlation in the left and the sub correlation in the right; and dotted-dashed 
lines show the result in the Cygnus region (Kohyama et al. 2010). The 41 GHz radio continuum emission 
map was not used by Hibi et al. (2006) or Kohyama et al. (2010) to separate the two correlations they 
identified. The contour map includes the same parameters as those described in Figure 2. 
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Table 5. Fitting results for the color-color correlations obtained by a power-law function (Equation 

Case Slow case Steep case 

easel a = 0.659 ±0.003 a =1.03 ±0.02 

weighted least squares fitting b = -0.767 ± 0.003 b = -0.508 ± 0.004 

(I P I > 30°) 

case II a = 0.24 ±0.02 a = 0.88 ±0.07 

un- weighted least squares fitting b = -1.36 ± 0.02 b = -0.627 ± 0.001 

(1/3 I >30°) 

Hibi ct al. (2006) a = 0.65 a = 0.93 

Galactic plane & = - . 78 & = - . 56 

(| b | < 5°) (main correlation) (sub correlation) 

Kohyama et al. (2010) a = 0.73 a = 1.09 

Cygnus region b = —0.75 b = —0.44 



(80° <Z<90°, -4° <6<8°) 



Errors in this study are shown in la. 



Table 


6. Calculation methods for Ay. 


Notation 


Method 


Ay (case I) 


our method, assuming case I. 


Ay (case II) 


our method, assuming case II. 


Ay (SFD) 


Schlegel, Finkbeiner, & Davis (1998) 


A v (2MASS) 


Dobashi (2011) 
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dust temperature (case I) 




dust temperature (case II) 




Fig. 4. Dust temperature distributions [K]. Top and bottom columns show results for cases I and II, 
respectively. 
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Av (case I) 




Fig. 5. Galactic extinction distributions at V band [mag]. The extinction is shown in a common logarithm 
scale. Top and bottom columns show the results for cases I and II, respectively. 
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Av (case I) [mag] 

Fig. 6. Comparison between Ay of this study and SFD. The systematic difference is 5% and 21%, for Ay 
(SFD) < 1 and Ay (SFD) > 1 mag, respectively. The systematic difference for Ay (SFD) > 1 mag may 
reflect the difference in the spatial resolution of the dust temperature. The wider dispersion for Ay (SFD) 
< 0.1 mag may reflect an SFD assumption (see text). The contour map includes the same parameters 
as those described in Figure 2. The bin size is set to 0.1 in logarithmic units. Contour lines are drawn 
according to the number of data points, including 0.5, 0.1N, 0.5N, N, 2N, and ION in the diagram, where 
N = 1525 indicates the average number of data points. The solid line indicates a slope of unity. 
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2500 r 




•1.0 -0.5 0.0 0.5 1.0 
Av (SFD) - Av (case I) [mag] 



Fig. 7. Histogram of the Ay difference between this study and SFD. Bin size is 0.01 mag. The mean, 
median, and standard deviation are —0.28, —0.001, and 2.29 mag, respectively. 



Av (SFD) - Av (case I) 




-0.1 0.1 



Fig. 8. Maps of the Ay difference between this study and SFD. The range is ± 0.1 mag. 
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Av uncertainty 




Fig. 9. Map of the Ay uncertainty of our method estimated by Equation (5) [mag]. The value is shown 
in a common logarithm scale. 




Fig. 10. Area in which the significant difference between this study and SFD is confirmed. Black shading 
represents areas in which the difference is larger than the Ay uncertainty of our method. The area occupies 
28% of the region for | /3 \ > 30°. 



25 




Av (case I) [mag] 



Fig. 11. Comparison of Ay between this study and Dobashi (2011). The Ay (2MASS) are saturated at 
Ay (2MASS) > 20 mag. The dispersion in the vertical direction at Ay (case I) < 1 mag may be caused 
by the low S/N of Ay (2MASS) in high galactic latitudes. The solid line indicates a slope of unity. LMC 
and SMC are masked. The contour map includes the same parameters as those described in Figure 6. 
The average number of data points contained in each bin is 732. 
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-1.0 -0.5 0.0 0.5 1.0 
Av (2MASS) - Av (case I) [mag] 



Fig. 12. Histogram of the Ay difference between this study and Dobashi (20f 1 ). Bin size is 0.01 mag. 
The mean, median, and standard deviation of the Ay difference arc —0.58, —0.09, 3.69 mag, respectively. 
LMC and SMC are masked. 



Av (2MASS) - Av (case I) 




Fig. 13. Map of the Ay difference between this study and Dobashi (2011). The range is ± 0.1 mag. LMC 
and SMC are masked. 
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Fig. 14. Area in which the significant difference between this study and Dobashi (2011) is confirmed. 
Black shading represents areas in which the difference is larger than the Ay uncertainty of our method. 
The area occupies 81% of the region for | /5 > 30°. LMC and SMC are masked. 



Table 7. Results of Ay comparisons among this study. SFD, and Dobashi (2011). 





Effective 
region 


Wavelength 


Spatial resolution 


la Ay difference 
between this study 


Upper limit 
for Ay 






H 


[arcmin] 


[mag] 


[mag] 


Ay Dust temperature 


This study 


| P | > 30° 


100 & 140 


5 5 




over 100 


SFD 


All-sky 


100 & 240 


5 42 


2.29 


over 100 


Dobashi (2011) 


All-sky 


1.25 & 1.65 


1 - 12 


3.69 


20 



The Ay difference is compared with the result for case I; i.e., Ay (case I). 
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